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We investigate how much amount of Greenberger-Horne-Zeilinger (GHZ) entanglement is required 
in order to prepare a given multipartite state by local operations and classical communication 
(LOCG).We present a LOGO procedure that asymptotically converts GHZ states into an arbitrary 
multipartite pure state, whose conversion rate is given by the multipartite discord of the state. This 
reveals that the GHZ-entanglement cost of preparing a pure state is not higher than the multipartite 
discord of the state. It also provides an operational interpretation of multipartite discord for pure 
states, namely, the consumption rate of GHZ entanglement in the devised procedure. 
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I. INTRODUCTION 

Quantum entanglement lies at the heart of notable fea¬ 
tures of quantum physics and the power of quantum in¬ 
formation processing Its quantification is an issue 

of great importance. A possible approach is based on the 
idea that entanglement is a resource shared between dis¬ 
tant parties that cannot be generated by local operations 
and classical communication (LOCC). More specifically, 
a scale of entanglement can be created by considering 
asymptotic conversion (i.e., conversion of m copies of a 
state into n copies of another state for large n and m) 
under LOCC. Entanglement is then quantified as the op¬ 
timal rate of the asymptotic conversion, m/n, from a 
resource state to the given state (or vice versa). 

Bennett et al. showed that, in the bipartite case, 
the asymptotic conversion between any pure states is 
reversible so when quantifying entanglement of bi¬ 
partite systems, any choice of a pure entangled state as 
a resource state results in an equivalent quantification. 
The singlet state is a natural choice for a resource state, 
and the corresponding conversion rate is called the entan¬ 
glement cost (and it is called the distillable entanglement 
for the reverse direction) im. The entanglement cost 
of preparing a bipartite pure state has been shown to 
equal the entropy of entanglement or, equivalently, the 
von Neumann entropy of one of the subsystems Q. 

Entanglement of multipartite states can be quantified 
in the same way, but different choices of the resource 
state may give rise to independent quantifications. For 
example, there exist asymptotically inequivalent multi¬ 
partite states such as the Greenberger-Horne-Zeilinger 
(GHZ) state and the W state. Hence, one may choose a 
resource state and correspondingly define a multipartite 
entanglement cost as the optimal conversion rate from 
the resource state to a given state iini- However, even 
achievable rates (i.e., upper bounds for the optimal rate) 
are not known for any resource state, except in the case 
where the singlets shared among the multiple parties are 
used as resource [T^. This can be attributed to the dif¬ 
ficulty of analyzing general LOCC [l^. The LOCC con¬ 


version between multipartite states has been little stud¬ 
ied especially in the asymptotic limit, while there have 
been studies on the stochastic LOCC conversion [IJ, [T^ 
and the conversion under asymptotically non-entangling 
operations mill in the asymptotic limit. 

Here, we consider the GHZ entanglement cost of 
preparing a quantum state of an arbitrary number of par¬ 
ties, where the fc-partite GHZ state is defined as 

|GHZ) = ^(10)®'=+ 11)®'=), 

V2 

with |0) and |1) denoting orthogonal basis states for local 
subsystems. We present a LOCC procedure, £d, for the 
asymptotic conversion from GHZ states to a multipartite 
pure state if, namely, 

GH •0®’^ 

for sufficiently large n. We find that the conversion rate 
of this procedure, Rd, is given by the multipartite discord 
of state ip, captured by the relative entropy of discord. 
This implies that the optimal rate is upper bounded by 
the relative entropy of discord. Our study further pro¬ 
vides an operational interpretation of multipartite dis¬ 
cord for pure states, namely, the consumption rate of 
GHZ entanglement in the devised procedure. In the mul¬ 
tipartite setting, quantum discord of pure states is dis¬ 
tinct from entanglement (although they are equal for bi¬ 
partite pure states), but its roles and meanings are not 
well understood compared to those of bipartite discord 
in several information tasks [isl - l^ . 

II. DEFINITIONS 

We consider a multipartite system P = {Pi,... ,Pk}, 
consisting of k subsystems of an arbitrary dimension. A 
/c-partite quantum state is called fully separable if it can 
be written in a form of J2jPjPi,j G) • • • 0 Pk,j- We de¬ 
note the set of separable states by S. We also denote 
an arbitrary orthonormal basis for the j-th subsystem by 
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{|a;j)}, where xj's are integers from 0 to dim(Pj) — 1. 
We can then construct a separable basis for the whole 
system as {|a;i,..., a:fc) = |xi) 0 • • • 0 |xfe)}. The ba¬ 
sis states for the whole system contain no nonclassical 
correlation, so we call their classical mixtures classically 
correlated states. To clarify, a fc-partite quantum state a 
is said to be classically correlated if it can be written as 

denote the set of classically correlated states by C. 

Apart from the quantification of entanglement based 
on the LOCC conversion, there is an axiomatic approach 
to measuring entanglement. In this approach, an entan¬ 
glement measure is given as a function of the density 
matrix that satisfies some desirable properties. One such 
measure is the relative entropy of entanglement defined 

asm 


= min5(p||cr) ( 1 ) 

ctGo 


where S'(p||cr) = —tr(plog 2 cr) — S{p) is the relative en¬ 
tropy of p to CT and S{p) = — tr(plog 2 p) is the von Neu¬ 
mann entropy of p. When dealing with many copies of 
a state, it is desirable to consider the regularized version 
of Er defined as [ 1 ^ 

f^(p)= hm ^-3^, 

Using the set of classically correlated states, C, the rel¬ 
ative entropy of discord | 2 fil-l^ and its regularized ver¬ 
sions are defined as 

Vr{p) =mmS{p\\<7), T>g(p) = lim \ (2) 

o-ec t->-oo t 


respectively. From the definitions, it is evident that 
^ Dr and ^ £r. As already mentioned, 
= Er for any bipartite pure state, but they are gen¬ 
erally not equal for fc-partite pure states with fc > 2 . 
It has been shown that one can always find a sepa¬ 
rable basis {|a;i,...,Xfe)} such that the complete set 
of projectors {Ilj} = , Xfc)(xi,..., Xfej} satisfies 

T^Rip) = 5'(p|| IljpIIj) [2^. We then obtain 




min 

{nu 


[- X! V;) log2 


( 3 ) 


for a pure state ip. 

To address the GHZ entanglement cost, we define the 
asymptotic preparation as follows. We say that a LOCC 
procedure, denoted by £, asymptotically prepare ip from 
GHZ states at rate R if F{[iP®'^,C{GHZ®'^)) 1 and 

min —?► i? as n —>■ oo, where A(p, a) = try^ p^ a pi is the 
fidelity between two quantum states p and a [ 2 ^. 


It is based on the asymptotic equipartition property 
(AEP) [s^ which will be briefly explained. Next, we 
present a LOCC procedure that prepares T from a cer¬ 
tain number of copies of a GHZ state, say m copies. Fi¬ 
nally, we prove that F{'^,ip^'^) 1 and m/n — T>^(ip) 

for n —?► oo, so they collectively verify that the LOCC 
procedure asymptotically prepare ip at rate 'D'^iip). 


A. Introduction of an approximate state 

Here we introduce a pure state 'I' that approximates 
ip®'^ (i.e., F{ip®'^, 4') —>■ 1 for n —>■ oo). For simplicity, we 
only consider tripartite systems, but the generalization to 
any fc-partite system is straightforward. 

We begin with summarizing the AEP [s^. Consider 
independent and identically distributed random variables 
..., Each of the variables has the same proba¬ 
bility distribution p{x) and the Shannon entropy F[. We 
define a typical set Ae to be a set of sequences \x'‘ = 
that satisfy ^ 

Then, the AEP states that the typical sequences {x’’} 
contain most of the probability, and the size of the typ¬ 
ical set, |.4e|, is about 2^. The AEP is summarized as 
follows. For any e > 0, 

Pr [x'' e Ae] > 1 — e, 

lAI > (l-e)2'(^-^), (4) 

\Ae\ 

for sufficiently large 1. 

Let us assume that we want to asymptotically pre¬ 
pare a tripartite state ip, where the dimensions of the 
three subsystems are I, rh and n. We then consider 
the state (p = ip®* prepared in a tripartite system P = 
{Pi,P 2 ,P 3 } where dim(Pi) = ([)*, dim(P 2 ) = {m)* and 
dim(P 3 ) = (h)‘. This state can be written in a separable 
basis {|a;i,..., Xk)} for system P as 

\(P) = \Pj)®*= ^ C{xi,X2,X3)\xi,X2,X3) (5) 

Xi,X2,X3 


with some coefficients C'(xi,X 2 ,X 3 ). We note that 
{|a::j)}’s can be any orthonormal bases for each subsys¬ 
tem, and we do not choose any particular separable ba¬ 
sis at this point. The coefficients \C{xi,X 2 ,X 3 )\^ can 
be considered a joint probability distribution of random 
variables Ai, A 2 , and A 3 , so we set p{xi,X 2 ,X 3 ) = 
|C'(a;i,a; 2 ,a; 3 )p and 

H = |C(a;i,a;2,a;3)plog2 \C{xi,X2,X3)\^. (6) 


III. MAIN RESULT 

We now provide our main result as follows. First, 
we introduce a pure state 4' that approximates ip®'^. 


We may consider independent and identically distributed 
random variables X^^\ ..., for j =1, 2, and 3. A 
typical set is correspondingly defined as 

Ae = {(4,4,4) : |fog2P(4:4>4) - H\^e}. (7) 
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Then, n = t ■ I copies of tjjj which is prepared in a sys¬ 
tem Q = {Qi, Q2, Qa} = is represented 

using the typical set Ae and its complement as 

I^)q” = 

{x\ , , 

" ( 8 ) 

+ C{x[,x‘2,x‘3)\x[,x'2,x'3). 

(.x\ ) 3 ^ 2 ’^ 3 ) 

The approximated state 'll is defined to consist of only 
the terms corresponding to the typical sequences, so 

\'^)Q = Nr^ C{x[,x’2,xl)\x[,x'2,xl) (9) 

(x[,x‘2,x‘^)£Ae 

where = E(x',4,x‘|C(2;'i, 4)|2 is a normal¬ 

ization factor. As the typical set contains most of the 
probability, we see that the fidelity between T and ■0®" 
approaches unity as n —^ 00 . In addition, it follows from 
the AEP that the number of terms in the expansion of 
T, |Ae|, is approximately 2^^. 

At this point, we note a previous method for bipartite 
asymptotic preparation Q in relation to our approach. 
In Ref. , one party, say the first party, prepares 0 q"q/ 
at their site and send the compressed Q^-part (the com¬ 
pression uses the AEP) to the second party using the 
quantum teleporation protocol. At the end of the proce¬ 
dure, the two parties share an approximated state of the 
same form as T. This method using the teleportation 
protocol can be applied only to multipartite states that 
have Schmidt decompositions such as q;| 000 )+/ 3 | 111 ) [ 31 |. 
On the other hand, our LOGO procedure directly con¬ 
verts singlets into the state It enables us to consider 
the approximated state T in multipartite setting, which 
can actually be created by LOCC as shown in the next 
subsection. 


B. Preparation of the approximate state 

We now present the LOCC procedure £d that con¬ 
verts log2 \ Ae\ copies of a GHZ state into 'I', so 
Cd{GHZ'^^°S2A^\) = ^f. We also only consider tri¬ 
partite systems, and the generalization to any fc-partite 
system is straightforward. For convenience of mathe¬ 
matical descriptions, we change the notation as follows. 
For x\, X2, and 0:3 such that {x{,X2,x’^) G Ae, consider 
the following sets: {{x\,X2,x^^)}, {a:^}, {x^}, and {ccg}, 
whose sizes are |Ae|, a, /3 and 7, respectively. We re¬ 
place {x\,X2,x^^) with y where y € {0,1,..., |Ae| — 1}. 
Similarly, we do the same for other sets as 

-^> {/ : / e 0, ...,a- I}, 

{ 2 ^ 2 } G 0,...,/?— I}, (10) 

{4} -)■ {h:hG0,...,j-l}. 


Because /, g, and h are completely determined by y, 
we denote them by /(y), y(y), and h{y). We can then 
rewrite the state d' of the system Q as 

^ l-Ael-l 

|^')Q=Af, " Y Ciy)\f{y),g{y),h{y))Q. (II) 

y=0 

The state 'I' can be obtained from m = |"log2 lAe]] 
copies of a GHZ state ([a:] is the smallest integer not less 
than x) by following four steps of LOCC. Assume that 
m copies of a GHZ state are prepared in an ancillary 
system Q' = {Q'i,Q2, Q3}, where dim(Qi) = dim(Q2) = 
dim(Q3) = \Af_\. We can rewrite GHZ®'^ as 

|GHZ)^r = 2 “^ Y^ \zi ■ ■ ■ Zm)\zi ■ ■ ■ Zm)\zi ■ ■ ■ Zm) 
2^-1 

= 2"^ Y |y:2/i2/)Q') 

y=o 

where y is the decimal representation of binary strings 
zi - ■ ■ Zm- Because \Ae\ < 2™, we may discard the terms 
other than those with 0 < y ^ \Ae_ \ — I (by a simple local 
operation), so we have 

^ l-Ad-l 

Y \y^y^y)(i'- 

y=o 

The first step is to change the coefficients from 2 “^ 
to Nr^^‘^G{y) as 

1.47-1 ^ 1 . 47-1 

Y l2/-y>y)Q'" Y Giy)\y,y,y)Q>. 

y =0 y^O 

This can be done by local operations of any party. Con¬ 
sider a local measurement described by the measurement 
operators 

_j|.A7-l I 

m,=n:^ Y ^(y©i)|y>(i/l, |A|-i[ 

M=0 J 

( 12 ) 

where 0 denotes addition modulo \Ae\- It is easy 
to check that they satisfy the completeness relation 
After the measurement, if the outcome 
of the measurement is j , each of the parties applies a uni¬ 
tary operation U: \y) —>■ |y 0 y) to complete the first step. 
The second step is to create the state \f{y),h{y),g{y))Q 
in the system Q, which is initially prepared in |000)q. 

^ I-A 7-1 

Y C'(y)|y,y,y)Q' 0 |/(y),/l(y),y(y))Q• 
!/=o 

This is achieved by a local unitary operation on Q\Q\ 
that transforms \y)Q'^ 0 |0)qi into \y)Q'^ 0 |/(y))Qi, and 
similar local unitary operations on Q2Q2 and QsQ's- 
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The third step is to disentangle Q 2 and Qg as 

C'(y)lj/)Q'i ® I/(J/)>5(2/),^(2 /))q- (13) 

v=o 

In Oder to perform this step, the second and third parties 
perform local measurements on their systems using the 
measurement operators 

{M, = \Ar^j\j){j\j^}, 

where J is a complex Hadamard operation, defined as 
{y\'^\y') = exp[27ri • yy' /\Ae\]- Depending on the mea¬ 
surement outcomes, the first party can perform a phase- 
shifting operation to complete the third step. 

The final step is to disentangle the ancillary system Q'^ 
from Q to obtain as 


phase-shifting operation by the first party that trans¬ 
forms 

\f{y)) -t exppTTi jP'if{y)]\f{y)). 

Similarly, consider phase-shifting operations by the 
second and the third parties that transform 

\g{y)) exp[27rt j-fg{y)]\g{y)), 

\h{y)) -)> exppTTi jh{y)]\h{y)). 

Applying those operations completes the fourth step. 


C. Asymptotic preparation from GHZ states at a 
rate equal to quantum discord 


^ lAI-i 

c'(y)l/(y)’ff(y)’^(y))Q = I^)q- 

v=0 

(14) 

To address the procedure for this step, we introduce a 
variable K(y), which is defined as 

K{y) = Pjfiy)+ -fg{y) + h{y). (15) 

Then, we can consider a linear isometry transformation 
V on the system Q[ which transforms \y) into \K{y)). 
This requires additional ancillary qubits in Qi so that 
the dimension of Q[ equals Ar(|Ac| — 1) -I- 1. Applying V 
to the state in Eq. (USD gives 


® \fiy)^ 9 iy)^Hy))Q- 


y=0 


Next, the first party performs a measurement on the an¬ 
cillary system Q[ using the operators 


|o)(jVt 


[k{\aa-i) + i ]2 


0< j ^ A(|A| -1) 


where J is a complex Hadamard operation defined as 


We denote hy Cd the LOCO procedure described in 
the previous subsection. We have shown that it prepares 
4- from j 

In practice, it can asymptotically prepare ip from GHZ 
states at rate 22^. This can be explained as follows. 
First, the AEP implies that the fidelity F('I','0®"') ap¬ 
proaches unity as n —?► 00, so asymptotically pre¬ 
pares 'i/)®". Second, the conversion rate, which is given 
by n“^|’log 2 lAel], can be reduced down to 'D^{'ip). Re¬ 
call that the typical set in Eq. (O depends on the sep¬ 
arable basis {a;i,X 2 ,X 3 }, which we have not specified 
yet. We now choose the separable basis so that it gives 
the minimum size of the typical set. Let us denote the 
minimum size by log 2 |A*| and denote the corresponding 
Shannon entropy by H*. It then follows from the AEP 
that log 2 |A*| ~ IH*. In addition, a comparison between 
Eqs. m and (| 6 |) leads to H — 'Ey. Putting these 
together, we have n“^[log 2 |A*]| = t“^l“^|’log 2 |A*|] « 
t~^T>fi{il!®*). Finally, from the definition of 22^, the con¬ 
version rate is approximately found to be 22 ^ {ip) for large 
t. We thus reach the following theorem. 

Theorem. The LOCC procedure Cd asymptotically 
prepares a multipartite pure state ip from GHZ states at 
rate 'D’^{ip). Namely, for any e > 0, d > 0, 


{y\J\y') = exp [27rz yy' [K{\Ae\ - 1 ) -f 1 ] 

for 0 ^ 2 /, y' ^ AdAel — 1) -I- 1. The completeness rela¬ 
tion, Mj Mj = I, can be checked from the orthogonal¬ 
ity of the Hadamard operation, jJ'l' = [2L(|Ae| —1)-|-1] 2. 
If the measurement outcome is j, the resulting state is 

c'(y)(2V^I^(2/)) ® l/(y)>5(y)>^(2/))Q 

y=0 

IWI-1 

=| 0 )Qi E ™ ^\fiy)^9iy)^h{y)}Q■ 

y=o 

Finally, local phase-shifting operations can remove the 
phase exp[—27ri jK{y)] by using Eq. (fT51) . Consider a 


F(/:D(G22Z®ri°S2lA:il)^^®n^ > 1 -e, 

\n-^\log2\A*l]-V^{ip)\ <S. ^ 

for sufficiently large n. 

proof — The LOCC procedure, £_d, works for any sep¬ 
arable basis {\xi,X 2 ,X 3 )'\ in Eq. (O. In a given separable 
basis, the Shannon entropy of the random variables Ai, 
A 2 , and Ag is given as 

H = |G(a;i,a;2,a;3)plog2 |C(a;i,a:2,a;3)|^, 

where |G(a;i,X 2 ,a; 3 )P = tr(|a;i,a; 2 ,a: 3 )(a;i,a; 2 ,a; 3 |' 0 ®‘). 
We note again that the typical set Ae depends on the 
separable basis and so does the Shannon entropy. We 
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can make H equal by choosing a suitable sep¬ 

arable basis. Using the expression of Vn in Eq. , we 
can write 

min I - y^tr{\xi,X2,X3){xi,X2,X3\ 

{\xx,X'^,X3){x-i,X2,X3\} L 

xV>®*) log2 tr{\xi,X2,Xi){x\,X2,Xi\-4>'^^')\^. 

Therefore, by choosing a separable basis that attains the 
minimum in the above equation, we have H = 

For clarity, we add a superscript * to the corresponding 
typical set and Shannon entropy as A* and H*. Then, it 
follows from the AEP that, for any e' > 0 , 

> (1 - e'), 


for sufficiently large /. The fidelity between i/;®" in 
Eq. (0 and 4 - = ) in Eq. © is given 

as 




E 

(xba:^,a:l,)eA 


Nf, 



In addition, applying the first inequality in Eq. dni) gives 

1 

N^, > (1 — e') 2. By choosing e' such that e' < 1 — (1 — e)^ 
for any e, we obtain 

F{Cd{GH > 1 -e 


for any e. 

As the LOCC procedure asymptotically prepares ■i/'"' 
from GHZr'°82l-^:il^ the conversion rate is given by 
log2 I A* I with n = t I (the ceiling function [ ] can 
be ignored for large n). Using the second inequality in 
Eq. (IT 71 ) and H* = V^{'4’®*'), it is straightforward to 
show that 


log2 

I 




< e' 


log2(l - e') 

I 


By dividing the above equation by t and choosing t' such 
that e' — log2(l — e') < S t (there always exists such 
e' > 0 for any I, t, S > 0 .), we have 

n t 

Finally, approaches 22 ^ (V') as t —>■ oo by 

the definition of , and it completes the proof of the 
theorem. 


IV. APPLICATIONS AND EXAMPLES 

In this section, we present upper and lower bounds for 
the GHZ entanglement cost (equivalently, the optimal 
rate) and examine them for several examples. We also 
compare the LOCC procedure Cd with another proce¬ 
dure that uses singlets as resource. Finally, we discuss 
the application of Cd to general mixed states. 


A. Upper and lower bounds for the GHZ 
entanglement cost 

We denote by the GHZ entanglement cost. The 
existence oi Cd implies that is an upper bound for 
Ec- In addition, one can easily show that the regular¬ 
ized version of any entanglement measure £ that is non¬ 
increasing under LOCC is a lower bound for Ac, pro¬ 
vided that £°°{GHZ) = 1 and it satisfies the continu¬ 
ity condition, i.e., if E{p^'^,a®'^) —>■ 1 as n —>■ oo, then 
i|£(/9®”) — f (cr”)| — >■ 0 as n —>• oo. The relative entropy 
of entanglement defined in Eq. o satisfies the continuity 
condition and £^{GHZ) = 1 that leads to 

For generalized GHZ states, y^| 000 ) + y/1 — p|lll), the 
two bounds coincide, so that our procedure is optimal 
and Ec = = -p\og2P - {1 - p) log2(l - p). In ad¬ 

dition, it is known that GHZ states can be distilled from 
the generalized GHZ states at the same rate [ 3 l|. 

The GHZ entanglement cost can be either greater or 
less than 1 as £^ and can be so. For instance, the 
GHZ entanglement cost of generalized GHZ states is less 
than or equal to one as we have shown. For the state 
(|000) -I- I -I- ll))/\/2 where |+) = (|0) + |l))/-\/ 2 , we have 
£^ = 1 from > £Tb-.c + ^a-.b (see Ref. [M IH 

for the inequality). In addition, one can easily see that 
a single copy of the GHZ state can be converted to a 
single copy of the state (|000) + | + ll ))/-\/2 by LOCC, 
so Ec = 1. The regularized discord is not known, 
but Eb = 1.5 so that Ec = I ^ E'^ ^ 1.5. In the 
case of the W state, |IU) = (jOOl) + |010) + |100))/-\/3, 
£^ = 2 log2 3 — 2 Ri 1.170 dsl and E^ = log2 3 r; 1.585, 
so that 1.170 < Ec < 1.585. 


B. Comparison with preparation of states from 
singlets 

For comparison, we consider a LOCC procedure that 
prepares a multipartite pure state tpQ^ from the sin¬ 
glets shared among Q [l^. In this procedure, a single 
party, say the first party, locally prepares the pure state 
^q" and uses quantum teleportation Q and data com¬ 
pression to distribute the state. For the tripartite 
case, it requires nS2 singlets between Q1-Q2 and nS^ 
singlets between Q1-Q3, where Si denotes the von Neu¬ 
mann entropy of the i-th subsystem. Therefore, it re¬ 
quires n{S2 + S3) singlets in total. Considering all the 
permutations of Qi, Q2 and Q3, the consumption rate of 
singlets is 

Rt = Si + S2 + S2 — maxis'!, ^2, ^a}. 

Since a singlet is obtainable from a single copy of a GHZ 
state, this procedure can also be achieved by consum¬ 
ing GHZ entanglement at rate Rt- We have no proof 
for E^ ^ Rt for general states. However, it has been 
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shown in Ref. [s^ that Vr ^ Rt (so ^ Rt) for 
a few kinds of three-qubit states including the gener¬ 
alized W states q;|001) -I- /3|010) -I- 7 IIOO), and gener¬ 
alized GHZ states a|000) -I- /3|111). For instance, for 
the W state, Vn rs 1.585 and Rt ~ 1.837. We also 
realize that for the fc-qubit generalized 

GHZ states, a|0®^) -|- /3|1®^), where we have used Rt = 
[ ^i] ~ maxis'!,..., Sfc}. In the case of the fc-qubit 

W state (|00---l)-t|00---10)-f----k|100---0))/Vfc, one 
can check that Dr ^ Rt for any fc. 

C. Application of Co to general mixed states 

Before presenting how to prepare general mixed states 
by using Ljj, we note that the relation Ec ^ Rr does 
not hold for general mixed states. This can be shown by 
taking a counter example of a bipartite state, 

(1 - 2p)|$+)($+| +p(|00)(00 | + | 11 )( 11 |) 

with 1$+) = ^1/2(|00) -I- |11)) and 0 ^ p < 1/2. 
The entanglement cost of the state is Ec = 772(1/2-1- 
!/p(l - p)) [13 where H2{x) = -a:log 2 a; - (1 - 

a;)log 2 (l — x), and the relative entropy of discord is 
= 1 — 772 (p). One can find that Ec > Rr for 
any 0 < p < 1/2. 

For a general mixed state a, our procedure can be 
applied to prepare the pure states ipi^s such that a = 
J2iPi4’i- Then, the mixed state a can be obtained by 
classically mixing them. This technique has already been 
used to generalize the entanglement cost of bipartite pure 
states to bipartite mixed states Q- In our case, the rate 
of the GHZ entanglement consumption is given by 

Rd{p) = lim (18) 

t^OO t 


where 

roicr) = inf : a = ^Pii/’*| . (19) 

We see that the rate Rjj vanishes for fully separable 
states. For bipartite systems, discord VR{'ijji) in Eq. (fl^ 
is equal to the entropy of entanglement, so Rd{p) is re¬ 
duced to the regularized version of the entanglement of 
formation Q. It is known that the regularized entangle¬ 
ment of formation is equal to the entanglement cost for 
bipartite states 0 - 


V. CONCLUSION 

We have suggested a LOCC procedure that asymptot¬ 
ically prepares an arbitrary pure state from GHZ states 
where the conversion rate is found to be the multipartite 
quantum discord. It reveals that the GHZ entanglement 
cost of preparing a multipartite pure state is not higher 
than a multipartite quantum discord of the state. Our 
work provides an operational interpretation of multipar¬ 
tite quantum discord in relation to a multipartite entan¬ 
glement cost. 
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